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Abstract 

The performance of linear receivers in the presence of co-channel interference in Rayleigh channels is a 
fundamental problem in wireless communications. Performance evaluation for these systems is well-known for 
receive arrays where the antennas are close enough to experience equal average SNRs from a source. In contrast, 
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is possible for the two-user scenario. In this paper, we derive closed form results for the probability density function 
(pdf) and cumulative distribution function (cdf) of the output signal to interference plus noise ratio (SINR) and 
signal to noise ratio (SNR) of minimum mean squared error (MMSE) and zero forcing (ZF) receivers in independent 
Rayleigh channels with arbitrary numbers of receive antennas. The results are verified by Monte Carlo simulations 
and high SNR approximations are also derived. The results enable further system analysis such as the evaluation 
of outage probability, bit error rate (BER) and capacity. 
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I. Introduction 

A macrodiversity multiple input multiple output (MIMO) system is considered in this paper to denote 
a system where both the transmit antennas and receive antennas are widely separated. As a result, the 
slow fading experienced on all links is different and each link has a different average signal to noise ratio 
(SNR). There is considerable interest in such systems from a variety of perspectives. They arise naturally 
in network MIMO [QQ|, ED, (3]| and in other types of base station (BS) collaboration 01. A simplified 
model, namely the classical circular Wyner model [0, is widely used in network MIMO systems, but is 
too restrictive to be useful in modern macrodiversity MIMO systems such as those proposed in 3GPP 
LTE- Advanced standards [7|. In the circular Wyner model, it is assumed that a given cell only experiences 
interference from two adjacent cells and this interference has a fixed level given by a particular fraction 
of the desired signal power. In contrast, the general model discussed in this paper does not make any such 
assumptions and assumes as many interferers as the system permits with arbitrary powers [4J. Macro- 
diversity can also occur as a result of collaborative MIMO ideas [5, p. 69]. In addition, the fundamental 
channel model, where each link has a different SNR, is closely related to the MIMO channel models in 
|[8l. Note that the multiple distributed transmit antennas could correspond to many single-antenna users, 
one MIMO transmitter with distributed antennas or variations of the two. 

The performance of linear receivers in such a macrodiversity system has been investigated via simulation 
||9l , iflOl . but very few analytical results appear to be available. Hence, in this paper we consider an 
analytical treatment of signal to interference plus noise (SINR) performance for two types of linear 
receivers: minimum mean squared error (MMSE) and zero forcing (ZF) receivers. We focus on the baseline 
case of a flat fading Rayleigh channel, where the links are all independent but have different SNRs due to 
the geographical separation. This subject has been well-studied in the micro-diversity case IfTTTl . Ifl2l . Ifl3ll , 
Ifl4l where there may be distributed sources, but the receive antennas are closely spaced. The performance 
metric of interest is the SINR/SNR distribution since this also leads to results for bit-error-rate (BER), 
symbol-error-rate (SER), outage probability and capacity, etc. 

The difficulty in analyzing macrodiversity systems is that there is no coherent methodology currently 
available to handle the type of channel matrices that occur. In independently distributed Rayleigh channels, 
basic results in statistics have been used to great effect IfTTTl . Ifl2l . In the presence of correlation, if a 
Kronecker correlation structure is assumed, there are also many results available ifTiSTl . ifTTl . lfT8l . |fT9l . 
These results tend to have their origins in multivariate statistics GUI and make heavy use of hypergeometric 
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function theory [1211 . Unfortunately, no such theory seems to be available for complex Gaussian matrices 
where every element has a different variance, i.e., the macrodiversity case. As a result, we focus on a case 
where progress is possible; the two-user scenario. For this scenario, the sources are two single antenna 
users or a single user with two widely separated transmit antennas. The sources communicate with an 
arbitrary number of base stations each with a single receive antenna or a single base station with an 
arbitrary number of widely distributed receive antennas. In this scenario, user 1 is detected with user 2 as 
the interferer and then vice-versa. A particular example of this scenario is also considered, where a three 
sector cluster in a network MIMO system communicates with two single antenna users. For the general 
two-user scenario, we are able to derive the exact closed form SINR/SNR distribution for both MMSE 
and ZF receivers. In addition to the exact SINR/SNR analysis, we also derive high SNR results for the 
SER of MMSE and ZF receivers. These results lead to a simple metric which relates system performance 
to the average link SNRs and therefore provides insight into the effect of these SNRs. In particular, we 
establish the following key observations and results. An exact analysis for the dual-user case is obtained. 
The methodology is not extendable to the multi-user case which suggests that approximations and/or 
bounds may be needed to handle more than two users. Simple, high SNR approximations are provided 
for the SER. These novel expressions provide a functional link between performance and the channel 
powers which is more accurate than previous measures such as orthogonality. The analysis shows that the 
performance of MMSE and ZF receivers becomes more different when the users have "parallel" channel 
powers. At low SINR, performance is enhanced by diversity (when the desired user has roughly equal 
channel powers at the receive antennas), whereas at high SINR performance is enhanced by a subset of 
receive antennas having a high desired power and a low interference power. 

The rest of the paper is laid out as follows. Section II describes the system model and receiver types. The 
main analysis is given in Sec. III. Sections IV and V give numerical results and conclusions, respectively. 

II. System Model and Receiver Types 

A. System Model 

Consider two single-antenna users communicating with tir distributed receive antennas in an indepen- 
dent flat Rayleigh fading environment. The C nRXl received vector is given by 



r = Hs + n, 



(1) 
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Fig. 1. A network MIMO system with a 3 sector cluster. 

where n is the C URXl additive-white-Gaussian-noise (AWGN) vector, s = (si,S2) T contains the two 
transmitted symbols from user 1 and user 2 and H is the C nRX2 channel matrix. The complex transmit 
vector, s, is normalized so that -E 1 { I -^l 1 2 } = -^{l^l 2 } = 1- The Gaussian noise vector, n ~ CM (0, a 2 I), 
has independent entries with E {\rii\ 2 } = a 2 , for i = 1, 2, . . . , ur. The channel matrix contains independent 
elements, H ik ~ CM(0, P ik ), where E{\H ik \ 2 } = P ik . 

A particular example of this scenario is shown in Fig[TJ where three BSs collaborate via a central backhaul 
processing unit (BPU) in the shaded three sector cluster to serve two single antenna users. In Fig. [Q it is 
clear that the geographical spread of users and receivers gives rise to a 3 x 2 channel matrix, H, where 
all the P ik values are different. 

B. Receiver Types 

At the receiver, the n R distributed antennas perform linear combining. Hence, the output of the combiner 
is f = W H r, where W is an ur X 2 weight matrix. The form of W and the resulting output SINR/SNR are 
well known for MMSE and ZF receivers. These results are summarized below. Without loss of generality, 
we assume that the index of the desired user is i = 1 and we denote the first column of W by Wy. In 
practice, user 1 is the desired user with user 2 the interferer, followed by user 2 being detected in the 
presence of interference from user 1. From [fTTTl . [fT3l . lfT5l the combining vector and output SINR of the 
MMSE receiver are given by 

w 1 =R~%, (2) 
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SINR = hfR'%, (3) 

where R = h 2 h 2 + a 1 ! and h\,h 2 denote columns 1 and 2 of H. From lfT3~l . 0~4| the combining matrix 
W = H (H H H) 1 and output SNR of the ZF receiver for n R > 2 are given by 

SNR = p — , (4) 



a 2 



(H H H) 
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where [B] n indicates the (1, l) th element of matrix B. 

III. System Analysis 

In this section, we derive the CDFs of the output SINR/SNR of MMSE and ZF receivers. First, we 
present some useful results as follows. 

A. ZF Analysis 

Let Z be the output SNR of a ZF receiver as given in ©. Z can be written as lfT5l 

Z = — 2 h['Mh { . (5) 



where M = I - h 2 (h 2 h 2 ^j h 2 . The characteristic function (cf) of Z is [|20l, J25 
Note that M and h\ are independent and the pdfs of hi and h 2 are given by 



(6) 



f{hk) = ^W 6 "^' 1 ^' for fc=i ' 2, (v) 

where the n R x n R matrix P k = diag(P lfe , P 2ik , . . . , P nR ,k)- Next, by using Lemma [2l the cf conditioned 
on h 2 becomes 



1 

\i-ji\MPA' 



^(^2) = - T ,,, /nr (8) 



Substituting M in © and rearranging gives 



>*(*l*2) = x , (9) 
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Fg(z) = cr 2 y^y^ fc Ji (nik,di,m ik ,fa,z) +ipikl2 (fiik,ai,rhik,fa,z) +^13 (n ik , &i, fh ik , fa, z^j . (10) 



i— 1 fc^z 



Z 2 (a, &, c, d, x) 



I\ (a, 6, c, d, x) 



dx 



be — ad 



, - be \ _ h / adx 

In [ — - e bx e ° Ei[ 

ad 1 \ c 



E 1 (bx) 



(ID 



c(bc-ad) (bc-ad)' 



e~ bx e k 



-E 1 



( adx 

\ c J (be- ad) 



d , . d + 6 (fee — ad) ( be 

wEi (bx) 5 — In — - 

(be -ad) 2 \ad 

1 - e~ bx 



a (be — ad) 



(12) 



h (a, b, c, d, x) 



c(bc-ad) (bc-ad)" 



e~ bx e^E 1 



( adx 



a , a , / bc\ 

2 Ei 6a; 5 In — 

(bc-adf (bc-ad) 2 \adj 



a(l-e- bx ) l- e 



— bx 



be (be — ad) bed 



(13) 



where D = I — -\jtPi. In Appendix El the full cf is obtained by averaging the conditional cf in © over 
h 2 . Then, the resulting cf is inverted to give the final result, the cdf in (flOl) . In (flOl) . the coefficients, (p ik , 
tpik and uiik are defined in (1731) and the arguments, ra^, a^, m ii; and are defined in (l65l) . The integrals, 
fi (.), /a (•) and J 3 (.) are given in (j74T)-(f76l 



5. MMSE Analysis 

The complete MMSE analysis can be found in [|26l . In this paper, we repeat a few of the initial steps 
which are also needed for the later high SNR results. Let Z be the output SINR of an MMSE receiver 
given by ©. The cf of Z is given bv ifTTTl. 11251 



it) = E{e' tz } = E^ th " R " h ^ 



(14) 



As in the ZF analysis, we first obtain the cf conditioned on h 2 . The conditional cf is given by 

Mt]h2) = \I-j l tR-^[ 
Substituting R in (fT5l) and rearranging gives 

<t>z(t\h 2 ) = — 

\D\ (a 2 + h 2 D~ l h 2 



(15) 



(16) 
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2 \ ^ \ ^ 



F z (z) = a z y y ]ipikh \ cr 2 h ik ,a t ,m, ik ,^r,z \ + ip ik I 2 I a 2 h ik , a i: ih ik , -|, z ) + u} ik I 3 I a 2 n ik ,ai,m ik , z 



(17) 



l\ (a, 6, c, d, x) = 



1 



6c — ad 



, a\ b ( b 



e^i ( - + bxj - e bx e - E x l 



(18) 



I 2 (a, b, c, d, a;) = 



1 + dx 



(be -ad) c(bc~ad) 



e~ bx e 



1 ' )a I (1 + dx) a\ cd + be — ad a (a 

Ei 5- e " Ei - 

c 7 c(6c-ad) 2 



(6c — ady 



1 - e" 6 * 
a (6c — ad) 



(19) 



h (a, 6, c, d, x) 



(be -ad) c (be - ad) 



-bx , 



■Ei 



l + dx)a\ ad 2 + abd-b 2 c 



+ 



d 2 (be -ad)' 



-Ei \ - + bx 



(be — ad)^ 



re-Ei 



d 

1 - e- bx 
d (be — ad) 



(20) 



'Results pertaining to MMSE receiver. 



Comparing ([9]) with ( |16l) we observe that, as expected, the MMSE results converge to the ZF results as 
a 2 0. Again, Lemma [Dean be used to average the conditional cf in (fT6l) to obtain the unconditional 
cf and in turn the cdf of the output SINR of an MMSE receiver as in ll26l . For the sake of completeness, 
we give the final results in (fT7l)-(l20l). The cdf is defined in (flTI) in terms of I\ (.), I2 (■) and ^3 (•) which 
are given in (fT8l) -(T20l). Finally, the necessary constants ipik, ipik and Uik are given by 



Pik = Pa R 2 (a 2 r]ik + Airiik + Cik ~ (nR - 2)P i3 fj i kj , 

Uik — —^i2^i\ Vik- 



(21a) 
(21b) 
(21c) 



IV. SER Approximations 

In this section, we derive high SNR approximations for the SER of MMSE and ZF receivers. This is 
motivated by the complexity of the exact analysis and the importance of finding a simple, functional link 
between performance and the average link SNRs. 



A. ZF Analysis 

The conditional cf in © is a ratio of quadratic forms in h 2 . Hence, (f> z (t) = E{4> z (t\h2)} is the mean 
of a ratio of quadratic forms which can be approximated by the Laplace approximation ||27ll as 

**(*)« ^ ,T 2 1 y (22) 

Note that the second equality in (O follows from the result, E{u"Qu} = Tr(Q), where u ~ CJ\f(0,I) 
and Q is a Hermitian matrix. Expanding the denominator of (1221) gives 

which follows since .D and P 2 are diagonal and \D\ is the product of the diagonal entries of D. As the 
SNR grows, a 2 — > and keeping only the dominant power of a 2 in (12~3l ) gives 

Tr(P 2 , a 2 ^- 1 ) 
~ iiA£^ — _Z (24) 

Vzl J |Px|Tr (P^Pj (-j*)"*" 1 

Tr P 

Defining i? (P 1; P 2 ) = t — — x — - gi yes a metric which encapsulates the effects of the power matrices 

iPilTrfPj P2 



Pi and P 2 . For many modulations, the SER can be evaluated as a single integral of the moment generating 
function of the SNR flU. The mgf of the SNR is M z (s) = (f> z (~js). As an example, for MPSK the 
SER is flU, [|2! 

* 1 rT 



, M z (—^-de, (25) 
7T J n V sin U , 

where g = sin 2 (tt/M) and T = Substituting (01) in (EU gives the approximation 

P s « (G„7j " . (26) 
In (1261) . the average SNR is 7 = ^2 , and the diversity gain and array gain are given by 

G d = n R -l, G a =(${P x ,P 2 )Z) 
The constant integral, X, is given by 

i = I / T f^V"" 1 *. (27) 

If.n \ / 
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1 ( T ( 2n R -2 \ (-l)""- 1 "^ 2 k ( 2n R -2 \ sin (2 (n R -k-l)T) \ 

^gna-l 1 2 2 ("«" 1 ) V TlR - 1 J 2 2 "«" 3 ^ 1 ' \ h J 2 (n R - k ~ l) J ' 

(28) 



Note that the simple representation in (|26l) shows the diversity order of n R — 1 and the effect of the link 
powers on array gain controlled by the metric $ (Pi,P2). The integral X can be solved in closed form 
and the final result is given in (|28l) . 

B. MMSE Analysis 

By using the Laplace approximation for the expectation of the conditional cf in (fT6l) , we obtain 

W |D| (a 2 + Tr(£)- 1 P 2 )) 
As the SNR grows, a 2 — >■ and keeping only the dominant power of a 2 in (|29l gives 

n.R-1 



« ; , V 7 \ -. (30) 

VZW |P X | (Tr (Pr^-j*) 
As in the ZF analysis, the SER for MPSK can be approximated by 

where the diversity gain and array gain are given by 

G d = n R - 1, G a = ($(P 1 ,P 2 )l(P 1 ,P 2 )y 1/{nR - 1) . (32) 

In (O, X(Pi,P 2 ) is given by 



1 [ T g^-V sm 2nR 6 
n Jo 9o + sin 2 



X(P 1 ,P 2 ) = -/ * ; . 2 „ dfl, (33) 



where g = g Tr (P^P 2 ). The integral, X(P 1; P 2 ), can be solved in closed form by expanding the ratio 
of sin 2 6 terms in (|33l) as a polynomial and integrating term by term to get the final result in (|34|) . We 
note that, as expected, the diversity order of n R — 1 is observed in both receiver types and the difference 
only appears in the array gains. 

The approximate, high SNR result for ZF in (1261) is particularly useful since it is simpler than the MMSE 
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X(P 1 ,P 2 ) = (-1) 



X 




(34) 



version in (|3T| ), and at high SNR the performance of the two schemes is similar anyway. Hence, (|26l ) 
acts as a useful approximation for both ZF and MMSE and provides a remarkably compact relationship 
between SER and the link powers, via the single function, •§ (Pi,P 2 ). 



In this section, we derive exact high SNR results for the SER of MMSE and ZF receivers. The work in 
this section does not employ the Laplace type approximation used in (1221) and (|29l ) and hence produces 
exact asymptotics at the expense of increased complexity. The mathematical details are given in brief to 
avoid unnecessary detail. 

A. ZF Analysis 

The conditional cf in © is a ratio of quadratic forms in h 2 . As the SNR grows, a 2 — >■ and keeping 
only the dominant power of a 2 in © gives 



where o (.) is the standard "little-o" notation and represents the fact that only the dominant power of a 2 
is used in the approximation and 



V. High SNR Analysis 




(35) 



Hence, the unconditional cf, when a 2 — > 0, becomes 




(36) 
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K 



V TiPi 



i=l 



12 



(t u p„i + t v p u i / 



In 



Pul 

Pu2 



-In 



5a 

P«2 



Pul 
Pu2 



Pvl 
Pv2 



(39) 



Following the same mgf based procedure to obtain the SER as in Sec. [TV] we arrive at the following 
expression 

-G e d 



P., = (G (( n) + o (~f- Ged ) , (38) 
where, the diversity and array gains are given by 

G ed = n R - 1, G ea = \J<oTj 
The constant, K , can be found using Lemma 1 to obtain the final expression as in (1391 where 



-l/{n R -l) 



P; 



T, = - (40) 

llk^i* hi* i2 — -T il-T fc2 

5. MMSE Analysis 

Consider the expectation of the conditional cf expression in ([Tot . As the SNR grows, cr 2 — > and 
keeping only the dominant power of a 2 in (fT6l) gives 

^z(t) = ^o(j*)(^) " +o(a 2 ^" 1 )), (41) 

where 

^„ - 1 e < - 

\Pi\ \h^P^h 2 -s 



f<^) =*e\ (42) 



Following the mgf based procedure in Sec. [TV] the SER at high SNR becomes 

Y^y-, , W (43, 

From (|43l) , the approximate SER can be written in terms of the diversity gain and array gain as 

Ps={GealT G<ld + 0(r G ° d ), (44) 



12 



X e (P 1 ,P 2 



n R 

X f V T,P 



,i=i 



ng 



l — f; {$,p - 1, ^) - r^p ( 



where 



G ed = n R -l, G ea = (l e (P 1 ,P 2 )) 



-l/(n R -l) 



Using Lemma 1, K (— s) can be given by 

*•<-> = -/ 5S 

which can be simplified to obtain 

#oH0 = 



TV • 2 /A ^ii? — 1 



P 



sin 2 fl 



d0. 



\Pi + e 1 p 2 p 1 + e 2 p 2 \ 



d9 2 , 



o 



oo / p nf 

( ^ '"'2 1 / 



Pi 2 6 2 + Pji 



riR 

£ 



PjlPj2 



P 2 #2 + Pi 



d0 2 , 



.1=1 " "° "/ \i=l 

where Tj is given in (l40l) . Equation (|47l) can be solved in closed form to give 



i=i 



where 



*. = £ 



Pi2^iPtlPfc2 + Pfc2^fcPjlPi2 
Pfclp2 — Pb2pil 



(45) 



(46) 



(47) 



(48) 



(49) 



Substituting s — gj sin 2 ^ in (|48l and then substituting P ( — m ^451) and integrating over 9 gives 
the result in (l50l) where 



P fm, a) 



n xsin 2 ^ 



sin 2m d0. 



(51) 



Clearly the exact asymptotics, especially for the MMSE case, are substantially more complex than the 
approximations in Sec. [IV] Also, the relationship between SER and the link powers is far more involved. 
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S1NR [dBJ 



Fig. 2. Analytical and simulated cdfs for the output SINR of an MMSE receiver for scenarios S1-S5 listed in Table 1 at p = 5 dB with 
? = 1. 




SINR [dB] 



Fig. 3. Analytical and simulated cdfs for the output SINR of an MMSE receiver for scenarios S6-S10 listed in Table 1 at p = 5 dB with 

? = 20. 



VI. Numerical Results 



In this section, we verify the analysis by Monte Carlo simulations using the network MIMO scenario 
in Fig. Q] [26]. We also consider some special cases of Pi and P 2 in order to investigate the effect of the 
macrodiversity powers on performance. For the two-user system in Fig. CD we consider the desired user 
to be user 1 and parameterize the system by three parameters. The average received signal to noise ratio 
is defined by p = Tr (Pi) /tircx 2 . The total signal to interference ratio is defined by q = Tr (Pi) /Tr (P 2 ). 
The spread of the signal power across the three antennas is assumed to follow an exponential profile, as 
in [fTTI . so that a range of possibilities can be covered with only one parameter. The exponential profile 
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SINR 

Fig. 4. Analytical and simulated cdfs for the output SINR of an MMSE receiver for scenarios (S1,S6) and (S3,S8) in Table 1. 



TABLE I 

Parameters for Figures[2]and[3] 





Decay Parameter 




Sc. No. 


Desired 


Interfering 




SI 


a = 0.2 


a = 0.2 




S2 


a = 0.2 


a = 1 




S3 


a = 0.2 


a = 5 




S4 


a = 1 


a = 1 




S5 


a = 1 


a = 0.2 




S6 


a = 0.2 


a = 0.2 


20 


S7 


a = 0.2 


a = 1 


20 


S8 


a = 0.2 


a = 5 


20 


S9 


a = 1 


a = 1 


20 


S10 


a = 1 


a = 0.2 


20 



is defined by 

P ik = K k (a)a i -\ (52) 

for receive antenna i, source k where 

K k {a) = Tx{P k ) / (1 + a + a 2 ), k = 1,2, (53) 

and a > is the parameter controlling the uniformity of the powers across the antennas. Note that as 
a — > the received power is dominant at the first antenna, as a becomes large (a>l) the third antenna 
is dominant and as a — > 1 there is an even spread, as in the standard microdiversity scenario. Although 




7 [dB] 



Fig. 5. SER and high SNR approximations for MMSE/ZF receivers using QPSK modulation in flat Rayleigh fading for four arbitrary drops. 




7 [dB] 



Fig. 6. SER of a ZF receiver using QPSK modulation in Rayleigh flat fading for three arbitrary drops. 

we consider microdiversity (S4 and S9), this is in the context of exploring the effect of different P matrix 
structures. Physically, it is not sensible to directly compare microdiversity with macrodiversity as they 
involve different system structures. In microdiversity, there may be multiple users communicating with a 
single array at a single BS. In macrodiversity, the users may be communicating with distributed antennas 
located at different BS sites which are back-hauled together to enable joint transmission/reception. In 
Figs. |2ll3] we show cdf results for the ten scenarios given in Table HI In Fig. [21 we see that SI has the 
worst cdf since the sharply decaying power profile is identical for both desired and interfering source. 
Hence, there is reduced diversity, as most of the signal strength is seen at one antenna, and there is strong 
interference at each antenna. Scenario S3 is best at high SINR since in this interference limited situation 
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Tr(PiP 2 ) 

Fig. 7. SER of a MMSE receiver for QPSK modulation in flat Rayleigh fading at 7 = 35 dB. 

(5 = 1) it is best to have at least one antenna where there is minimal interference. This occurs with S3 
as the power profiles are opposing and the strongest desired signal aligns with the weakest interferes At 
low SINR, scenarios S4 and S5 are slightly better than S3 as they have full diversity (equal power at each 
antenna) which is beneficial in this SINR region. In Fig. similar results are observed with diversity 
being important at lower SINR (where S9 and S 10 are the best) and interference reduction being important 
at high SINR (where S8 is best). In Fig. 51 we consider the effect of q on two different macrodiversity 
scenarios in Table U In particular, we have shown results for SI and S6 and S3 and S8. As expected, SI 
and S3 have lower SINRs than S6 and S9 due to increased interference. However, SI is far more sensitive 
to ^ than S3. This is because S3 has opposing power profiles for the desired and interfering users so that 
the two sources are more orthogonal and interference plays a smaller part in performance. 
Next, we consider the high SNR results in Sees. [IV] and [V] In Fig. [5] the MMSE/ZF receivers are 
considered for four drops (Dl, D2, D3 and D4) of two users in the shaded coverage area of Fig. [Q Each 
user is dropped at a different random location (uniformly generated over the coverage area) and random 
lognormal shadow fading and path loss is considered where osf = 8dB (standard deviation of shadow 
fading) and 7 = 3.5 (path loss exponent). Hence, each user has a different distance and shadow fade to 
each BS and each drop results in a new P matrix. The transmit power of the sources is scaled so that 
all locations in the coverage area have a maximum received SNR greater than 3dB, at least 95% of the 
time. The maximum SNR is taken over the 3 BSs. For all four drops, the high SNR approximations from 
Sec. [IV] are shown to be very accurate for SERs below 1CT 2 , although for drop D2 the results are less 
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tight. Note that drop D2 has the greatest difference between the MMSE and ZF results. In general, the 
gap between MMSE and ZF results can be assessed by a comparison of (1271) with (1331) . Here, it can be 
seen that the approximate asymptotics are the same for ZF and MMSE when g — 0. Hence, scenarios 
where g is large, i.e., Tr (P^ 1 P 2 ) ~ 0, will create substantial differences between the two receivers. In 
Fig. |5J drop D2 had the smallest value of Tr (P^~ 1 P 2 ) and hence showed the greatest difference. Note 
that for Tr (P^~ 1 P 2 ) to be small, Pn ^> Pi2 is required for i = 1,2,..., tir. Hence, the power profiles for 
users 1 and 2 must be "parallel" in some sense, with any large value of P i2 aligning with an even larger 
value of Pa. In these "parallel" scenarios, MMSE and ZF results can exhibit greater differences. 
This methodology is used again in Fig. |6] for three drops and ZF results are shown. Again, the asymptotic 
results show good agreement at SERs below 10~ 2 . Furthermore, the difference between the approximations 
in Sec. |IV] and the exact asymptotics in Sec. |V] is shown to be minor. Hence the simple SER forms in ((261) 
and (|3Tb are particularly useful. Note that the power matrices, Pi and P 2 , are completely general, with 
the sole constraint being P ik > 0, Vi, k. As a result, it is likely that some combinations of powers can be 
found that will cause the approximate SERs to lose accuracy. However, in all the scenarios considered 
and all random drops simulated (see also [1261 ) the approximations have shown similar accuracy to the 
results in Figs. [5J and |6l 

In multiuser systems it is well-known that two users may be successfully detected if their channels are 
approximately orthogonal. In the context of a dual user system, where only the channel powers are 
considered, the analog would be that Tr(PiP 2 ) is small. To investigate this relationship, we generate a 
large number of random power matrices with a fixed total power. For user 1, the powers are generated 
by (1521 with a = 0.2. For user 2, the powers are independent uniform random variables which are scaled 
so that c = 1. For each pair, (Pi,P 2 ), we compute Tr(PiP 2 ) and the approximate SER of user 1 using 
using (T3T]). The results are plotted in Fig. |7J As expected, SER increases with Tr (PxP 2 ), although there 
is wide variation in the band of SER results. In comparison, the new metric, (Pi,P 2 ), has a one-to-one 
relationship with the approximate SER and carries far more information than ad-hoc measures such as 
Tr (PxP 2 ). 

VII. Conclusion 

In this paper, we derived the exact cdf of the output SINR/SNR for MMSE and ZF receivers in 
the presence of a single interfering user with an arbitrary number of receive antennas. To the best of our 
knowledge, this represents the first exact analysis of linear combining in macrodiversity systems. Although 
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not shown for reasons of space, the slightly simpler problem of obtaining the associated pdf can also be 
handled since the pdf expression in (|58l) is inherently computed en-route to the cdf in (|59l) . Numerical 
examples demonstrate the validity of the analysis across arbitrary drops and channels. This suggests that 
the analysis is also numerically robust. A high SNR analysis reveals simple SER results for both MMSE 
and ZF which provide insights into both diversity and array gain. It also provides a functional link between 
the performance of the macrodiversity system and the link SNRs. 

Appendix A 
Background Results 

The following lemma gives a compact method to calculate the expected value of a ratio of random 
variables with arbitrary integer powers |[23l . 

Lemma 1. Let U\, U2 and Z be three continuous random variables such that P(U2 > 0) = 1, Z > and 
Z = jjjt. Assuming that there exists a joint moment generating function (mgf) for U\ and U2, denoted 
■M-{6i, 62) = E (e 6,lf/l+e ' 2C72 ), then, for all positive integer values of m and n, the expected value of Z is 
given by 

1 r°° . F) m KA(0 -t\ 

dz. 



E > u ?\- 1 r z n-i dmM 



J (n-1)! J ' d6 r 



The following n-dimensional complex Gaussian integral identity will be used extensively. 

Lemma 2. Let A be an arbitrary n x n complex Hermitian positive definite matrix. Then, the following 
integral identity holds. 

/OO PCO 
... j ( x " Ax d.v x . ..(I.r, 
-00 J —00 



7T" 

LAI' 



where the complex n x 1 vector x = [x\, . . . , x n ] T , dxi = dxudxiQ, xu = Re(a;j) and xiq = Im(xj). 
Finally, we give the following partial fraction expansion from elementary algebra, 

- = V A (54) 

where j 2 = —1 and 

j\. = (55) 
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Appendix B 
ZF Analysis 

From E51 . the full cf can be obtained by averaging the conditional cf in © over h 2 . Using Lemma [TJ 
the full cf is given by 



1 



dE{( 



-01Z1-02Z2 



d9 1 



d6o 



(56) 



6»i=0 



where z\ = h 2 h 2 and z 2 = b^D x h 2 . Note that the dummy variables are 9i, 6 2 > and we have used a 
slight variation of Lemma [TJ in which the joint mgf has negative coefficients, for convenience. Using the 
pdf of h 2 in © to evaluate the expectation in (|56l) and using the Gaussian integral identity in Lemma |2] 
and a few simplifications, we obtain the following result. 



f°° d 



\D 1 -jtD 2 \ 



d9 2 , 



(57) 



6»i=0 



where D 1 = I + 6 l P 2 + 6 2 P 2 and D 2 = j i -P 1 (J + ^P 2 ). From the pdf and cdf of Z are 

i r°° 

fz(z) = 7T / <Pz(t'K Jtz dt, (58) 

and 

i pz /*oo 

F Z {z) = J <i> z (t)e-> tx dtdx. 



(59) 



By substituting (1571) into (1581) and (|59l) multiple integral forms for the pdf and cdf of Z are obtained as 



1 

2^ 



OO /"OO 



5_ 

^ 



d9 2 dt, 



(60) 



6»i=0 



^0*) = -nZ 

Z7T 



oo 3 



-jtx 



\D X - jtD z 



dOodtdx. 



(61) 



0i=O 



'0 J-ooJQ 

Since D\ and Z) 2 are diagonal, we can further simplify the expression in (|6"TT) with the substitutions, 
di = 1 + ^iPi2 + #2-P«2 and 6j = ^-Pii (1 + 0\Pi2)- Then, the integrand in (|6TT) , before differentiation, can 
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dJi 



8i=0 



i=l 



(n R - 2)P l2 ) 



p n R- 



hrhik) 



(64) 



Z POO POO 



2k Jo J- 



00 JO 



d6 2 dtdx. 



(62) 



Since the limits of integration are independent of 6\, we can interchange the order of differentiation and 
first perform the integration over t. To obtain this integral, we use the partial fraction expansion of J 
from (|54l) and apply the following integral identity from E4l . 



-jpx 



09 - j*)" 



dx 



2wp v - 1 e- 

r» 



p > 

P < 0, 



> 0, Re(P) > 0] 



This gives the result, 



jo 



fflx 



d6 2 dx, 



(63) 



where J\ = YH=i^i e Differentiating J\ term by term and setting 9\ 

(Pa. - Pki) {Pa + Pk2) and 5 ik = (P a - P kl ) P l2 P k2 . Also 

flik = (Pit — Pkl) , mik = {PilPk2 — PklPiV 



gives (1641) where 7. 



ik 



a 



a 2 P i2 



&i p > Pi — p ■ 
Pi Pi 

Using (1541 ), the product in the denominator of (|64l ) can be expanded as 

1 _ \ ^ Vik 



(65a) 
(65b) 



(66) 



where 



m 



Vik 



ik 



Yli£i tk (narhi k - n ik m a ) ' 
Substituting (l66l) in (|64l) gives (1681) where the constants are given by 



(67) 
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dJi 

8b\ 



R 



1=0 



= 1 k^i 



(n lk + 9 2 m ik ) (n ik + 9 2 m ik ) (n ik + 9 2 m ik f 



(68) 



2=1 fc 7^/ 



x6 2 <T 2 P? 2 Pl R -% k 



Ml '/ifcS 



7 A: 



(n ik + 9 2 m ik ) (n ik + 9 2 m, ik ) {n ik + 9 2 m ik y 



d9 2 dx. (72) 



^ _ (PjlPj2Pk2 — PklPk2Pj2) 

{Pi\Pk2 — -P«2-Pfcl) 

~ (Pn - P k l) (PilPi 2 - PklPa) nn \ 



{Pi\Pk2 — Pi2Pkl) 



Qk = m ik } — — ; - (71) 

Substituting from (|68l ) in to the cdf expression in (|63l) gives (T721 ). The desired cdf in (1721) is 

rewritten in (flOl) where 

<^ifc = ^™ fl ~ 2 (A^jfc - (n fl - 2)P i2 fj ik + Ofc) , (73a) 
Ak = r}dikPii R '\ (73b) 
= -Pa R ~%kP- 2 (73c) 

Note that when n R = 2, Q k = and 77^ = 1 for all i, k. The cdf in (flOl ) contains three types of double 
integral defined by 

r-x r-oo ^—bt-dtQ 

h (a, b, c,d,x) = / -d6dt, (74) 

Jo Jo a + c9 

r-x r-oo ^—bt—dte 

I2 (a, 6, c,d,x) = I / ^d9dt, (75) 

Jo Jo (a + c6) 

r-x r-oo ±r\ —bt—dt8 

I 3 (a,b,c,d,x)= / ——ffldt. (76) 

Jo Jo a + c9 

Each double integral can be evaluated using standard methods in terms of a sum of exponential integral 
functions as shown in (fTTI)-(fT3~l). where E\{x) = ^-dt. An outline of the solutions of £74])-(|76]) is 
given as follows. I\ in (|74|) can be solved by integrating over t first and making use of the two identities 
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in B21 



Aa 



oo -xe f /\ a \ 

—e d0 - V Bi (t) • <77) 

cfe In {p/ a) 



(78) 



o (x + a) (x + 0) f3 — a 

to solve the resulting integral over 9. Next, we note that I 2 follows directly from Ii as I 2 = — ^p. Hence, 
by differentiating the expression for Ii in (fTTT) . we obtain (fT2l) . In order to differentiate the exponential 
integral, we use Leibnitz's integration formula to give 

9 [Ex (aa)} = _ (79) 

da a 

I3 in (|76l) can also be solved by a similar approach as in ii and making use of the following integral 
identity from ll24l where necessary: 

-bt 

'0 
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